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Abstract 



We make use of continuum elasticity theory to investigate the collective 
modes that propagate along the edge of a two-dimensional electron liquid or 
crystal in a magnetic field. An exact solution of the equations of motion 
is obtained with the following simplifying assumptions: (i) The system is 
macroscopically homogeneous and isotropic in the half-plane delimited by the 
edge (ii) The electron-electron interaction is of finite range due to screening by 
external electrodes (iii) The system is nearly incompressible. At sufficiently 
small wave vector q we find a universal dispersion curve u> ~ q independent of 
the shear modulus. At larger wave vectors the dispersion can change its form 
in a manner dependent on the comparison of various length scales. We obtain 
analytical formulas for the dispersion and damping of the modes in various 
physical regimes. 
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I. INTRODUCTION 



The dynamical behavior of the edge of a two dimensional electron gas (2DEG) in a 
strong magnetic field is receiving considerable attention. The interest has been spurred by 
recent experiments in which an unusual tunneling current- voltage relation of the form 
/ = V 1 ^ , where V is the bias voltage and vq is the bulk filling factor, has been observed. 
Such a distinctly non-Fermi-liquid-like behavior can be explained @-f| on the assumption 
that there is only one branch of collective edge wave excitations (the so-called charged mode), 
and that the tunneling electron is initially accommodated as a coherent superposition of such 
waves spreading in time according to hydrodynamic equations of motion. 

Unfortunately, a complete theory of edge excitations is not yet available. For example, 
there is considerable controversy about the existence of additional "neutral modes", and 
whether they contribute to the tunneling characteristics or not |jf-[J . 

Even the more conventional charged mode has not been fully analyzed yet. The two most 
successful theories so far PJT0|| focus exclusively on the long-range part of the Coulomb 
interaction and should therefore be viewed as macroscopic versions of the random phase 
approximation (RPA). A complete description of the long wave dynamics requires, however, 
the inclusion of short range forces that arise from the relative motion of adjacent parts of the 
fluid. These can be pressure forces arising from the bulk modulus of the quantum mechanical 
fluid, as well as shear forces arising from more subtle positional correlations. Shear forces 
are, of course, essential in the crystalline phase, but they are also present in the liquid phase 
at non zero frequency ||11|| . In addition, there are viscous forces - completely ignored in the 
RPA - that cause damping of the collective modes. 

These post-RPA effects are well established in the theory of the dynamics of the uniform 
electron gas, where they are usually described in terms of "local field corrections" (see [|IT|j ). 
For purely longitudinal bulk modes these corrections are not too important as the direct 
Coulomb interaction controls the physics of long wavelength longitudinal fluctuations. The 
situation can be quite different for edge modes, or even for bulk modes in a magnetic field, 
because in these cases longitudinal and transverse channels are strongly coupled, and one 
cannot talk of purely longitudinal or transverse modes. Indeed, we have recently shown 
@] that the collective dynamics of a uniform electron gas, in the limit of infinite magnetic 
field, must necessarily include a shear force term, otherwise the frequencies of the collective 
modes vanish. This led us to the description of the electron g clS clS 8b continuum elastic 
medium, characterized by elastic constants (the bulk and the shear modulus) which control 
the dynamical response, and by viscosity coefficients, which control dissipation. 

Encouraged by the qualitative success of the continuum elasticity approach to the dy- 
namics of the uniform electron gas, in this paper we present its application to the problem 
of edge dynamics in a magnetic field. 

On a macroscopic scale, our system is modeled as a uniform continuum in the half-plane 
delimited by the edge. The external potential, that confines the electrons to the half-plane, 
is constant in the bulk of the system, and rises sharply at the edge. A great mathematical 
simplification follows from the assumption that the electron-electron interaction is of finite 
( microscopic) range. This assumption is justified for structures in which the long-range 
Coulomb interaction is screened by metallic electrodes. The short range of the interaction 
results in equations of motion for the elastic displacement field that are linear differential 
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equations, as opposed to integral equations in the approach of |],|T(| . 

The distinctive feature of our theory is that the presence of the edge is taken into ac- 
count via boundary conditions. In an infinite system, the solutions of the equations of motion 
are purely longitudinal or transverse plane waves characterized by a real wave vector. The 
presence of the edge allows the existence of solutions with an imaginary wave vector perpen- 
dicular to the edge, which therefore vanish exponentially as one moves into the bulk. It is 
evident that these "bound" solutions exist independently of the nature of the bulk, in par- 
ticular, irrespective of whether the bulk is incompressible or not. In addition, these solutions 
must have the property that the elastic stress vanish at the edge. Because the boundary con- 
ditions are linear, the problem can be solved exactly by elementary techniques, and analytic 
results for the dispersion and damping of the collective modes can be obtained. (The quali- 
tative changes brought about by the long range of the Coulomb interaction will be discussed 
in the Appendix). A particularly elegant solution can be obtained for nearly incompressible 
systems, in which a large electrostatic charging energy strongly opposes density fluctuations. 

Our main result is that, at sufficiently long wavelengths and high magnetic field there 
is only one low-frequency edge mode, whose dispersion depends on the strengths of the 
confining electric and magnetic fields, but not on the bulk and shear moduli of the system. 
The nature of this solution does not change when the system has crystalline order. At shorter 
wavelengths, the dispersion can change its form in a manner dependent on the comparison of 
various length scales. In particular, we point out the possibility of a crossover from a linear 
dispersion, controlled by the confining electric field, to a quadratic dispersion, controlled by 
the shear modulus of the electronic system, if the latter is sufficiently large. 

Our results support the idea that the dynamics of a sharp edge in a 2DEG is completely 
dominated by a single charged mode, with any additional structure associated with peculiar 
characteristics of the system (such as the fractional quantum Hall effect) becoming irrelevant 
at long wavelengths. 



II. THE MODEL 

We consider a two dimensional electron fluid on a uniform background of positive charges, 
in a constant magnetic field B = Bz. The system extends indefinitely for x < while 
presenting a sharp edge at a; ~ 0, parallel to the y axis. Since the width a of the edge 
is much smaller than the characteristic wavelength of the modes under consideration, we 
model the edge setting a ~ 0. The mass density, at equilibrium, is a constant p for x < 0, 
and vanishes for x > 0. Consistent with this assumption, the external electrostatic potential 
V ex t(x), created by the background of positive charges and by nearby gates, is taken to be 
constant for x < and to rise sharply at x = with a derivative 

_ e ( dV ext \ , . 

9 ~ m { dx ) x=0 

Here — e is the effective electron charge, i.e., the bare charge divided by the square root 
of the static dielectric constant of the host semiconductor, and m is the effective mass of 
the electron in the host semiconductor. The quantity g is the acceleration imparted by the 
external potential to an electron at the edge. The effect of this acceleration on the collective 
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edge modes is analogous to the effect of gravity on surface waves in a liquid - hence the 
notation. 

To describe the dynamics of the system we introduce the displacement field u(x, y, t) of 
the infinitesimal volume element at point (x, y) from its equilibrium position. This obeys 
the linearized equation of motion of continuum elasticity theory ||12|| : 



_ = _£a,^+ We (ax-), (2) 

where i,j are Cartesian indices, di is the derivative with respect to r iy u c = eB/mc is the 
cyclotron frequency, and 

Gij = KV ■ u5ij + fi(diUj + djUi - V • uSij) (3) 

is the elastic stress tensor with bulk modulus K and shear modulus /i. The first term on 
the right hand side of eq. (^) is the force exerted on the volume element by the surrounding 
medium, the second term is the Lorenz force. In writing these equations, we have assumed 
that the electron-electron interaction v(r — r') is of finite range in space. This means that 
v (q = 0) = / v(r)dr is finite and is included as part of the bulk modulus: 

K = v(q = 0)n 2 + K, (4) 

where K is the "proper" contribution arising from the kinetic and exchange-correlation 
energy and n = po/m is the particle density. This procedure is, of course, only justified 
when the spatial variation of the density is small over distances of the order of the range of 
the interaction d, so that qd -C 1, where q is the wave vector. We must also have qa -C 1 for 
the "sharp edge" description to be valid. Clearly, both conditions are satisfied at sufficiently 
long wavelengths. 

Substituting eq. ([3]) into eq. (0) we obtain the standard form of the equation of motion 



d 2 u 

dt 2 



C 2 V(V • u) + CfV 2 u + co c z x u, (5) 



where 



C 2 t = ± (6) 
Po 

is the square of the transverse sound velocity in the absence of the magnetic field and 

C? = - = C*-C$ (7) 
Po 

where Cf is the square of the longitudinal sound velocity in the absence of a magnetic field. 
As stated in the introduction, we assume that the system is nearly incompressible, in the 
sense that the charging energy v(q = 0)n gives the dominant contribution to the bulk 
modulus, and is much larger than any other energy scale in the problem, such as p/n or 
hu c . This implies that C\ ^> C t so that the difference between C 2 and Cf can and will 
be ignored in the following. We emphasize that this assumption (which is expected to be 
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reasonable for an electrically charged system) simplifies the calculations, but is not otherwise 
essential. 

To complete the definition of the model we must now specify the boundary conditions 
on the solution of eq. (H). In elasticity theory, the normal form of the boundary condition 
at a free surface is ( |p~2f ) 

E^i| ed9e % = 0> (8) 

where i = x or y, n (= x in this case) is the unit vector perpendicular to the edge, and the 
subscript "edge" means that the quantity on the left hand side must be evaluated at the 
position of the moving edge, that is, at the point of coordinates (u x (0, y), y). The physical 
significance of this condition is that there is no matter beyond the edge to exert a force on 
the system. 

This boundary condition must be slightly modified here because the external potential 
produces a stress in the edge region even when the system is in equilibrium. Let us denote 
by af^ this equilibrium stress. Then the free boundary conditions take the form 

Yl[<7i?(u x ,y) + (7^(^,2/)]% = 0. (9) 
j 

From the symmetry of the problem one sees that only the xx component of the 
equilibrium stress tensor is non zero, and from the equilibrium condition da x °J(x)/dx + 
e(po/m)dV ex t(x)/dx = one sees that 



'dai°) 



/ x=0 



gpo, (io) 



with g defined in eq. ([!]). Finally one is free to choose o~ x °J(x = 0, y) = so as to satisfy the 
free boundary conditions at equilibrium. 



Expanding eq. (0) to first order in u, and making use of eq. fllPf ), we finally obtain explicit 
forms for the two boundary conditions: 



2 , , 2 du x (0,y) du y (0,y) 
Ci V ■ u(0, y) + C t ( — 7^ ^- — ) +gu x {0,y) = (11) 

du y (0,y) | du x (Q,y) _ Q 
dx dy 

Together with the condition that the displacement field vanish for x — ► — oo, eqs. ([5]) and 



(p|) - (|12l) completely define the mathematical problem under consideration. 



III. EDGE WAVES IN A MAGNETIC FIELD 

Ignoring at first the boundary conditions it is easy to see that the solutions of eq. @ 
that vanish for x — > — oo can be chosen to have the form 

u(x,y,t) = ue iqy+Xx - iuJt , (13) 
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where q is a real wave vector parallel to the edge and A is, in general, a complex number, 
whose real part must be positive in order to ensure decay in the interior of the system. 

Taking the divergence and the curl of the equations of motion (|5|) and making use of 
eq. dl3|) , we obtain: 



2/ \ 2 2 
_, v _ ^ 

2 /-(2/\2 



-u; 2 - Cf (A - g 2 )] V • u - ioouj c (V x u) 



-a; 



^-?1(Vxu) 




0. 



(14) 
(15) 



These two linear homogeneous equations are compatible if and only if the square of the 
wavevector A 2 has one of the two values 



\ 2 - a 2 - ' ^ l ' 4- 



\ 



'u 2 {Ci + Cf)Y uj 2 (uj 2 -uj 2 ) 



2C?C? 



C?C? 



(16) 



In the limit C\ C t the two solutions for A 2 (which we denote A 2 and A 2 ) take the simple 
form 



X 2 - n 2 
A l — Q 



LU 2 - UJ 2 



A 



Q 



uj- 



(17) 
(18) 



With a little algebra, it is possible to calculate the corresponding eigenf unctions in the 
limit C\ — > oo: 



u; oc [x + 2sgri(g)y]e 



iqy+\q\x—iu}t 



(19) 



u t oc [x + isgn(q)yfl - UJ 2 /Cfq 2 y]e iqy+lq ^ 1 ~ uj2/c ' q2x - iujt . (20) 

Let us now turn to the problem of satisfying the boundary conditions (fP2| ) at x = 0. 
This can be done by forming a suitable superposition of the two independent solutions Ui 
and u t , namely 

u = aui + 6ut (21) 

where a and b are complex coefficients. In order to implement the boundary conditions of 
eq. (0) we must calculate the limit for C\ — > oo of the products CfV • uj(<). To accomplish 
this, we need to refine our calculation of the eigenf unctions u^ t ) by including terms of order 
1/Cf. This can be done straightforwardly, with the help of eqs. ([T7|) and flT8|), and the 
result is 

limc^V ■ 0,(0, y) = ^-^"^ e .,y ( 22) 

Q 

and 

Umc^ooCfV ■ u t (0, y) = . (23) 
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Substituting eqs. fl22"|) and (|23|) in the boundary conditions (|ll])-(|r|) yields the following 
set of linear homogeneous equation for the coefficients a and b 

uuj c - u 2 sgn(q) 2 



[<? + — + 2C t 2 |g|]a + [ 9 + ^£ + 2Clsj\ - u>/C?q*\q\]b = 

2a+[2-uj 2 /C 2 q 2 }b = 0. (24) 

These two equations are compatible if and only if the frequency u satisfies the algebraic 
equation 



where 

X = 1^2 (26) 



\q\Q 

UJr 



z = w, (27) 
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« = (28) 



The complete solution for the displacement field is 



u(x, y, t) = { [(2 - £ 2 )e^ - 2e'«l V 1 ^]* 



+is#n(g)[(2 - £ 2 )e k|x - 2^1 - ^Mv^^y}^"^). (29) 



Equations ( p5|) and (|29f) are the central results of this paper. In particular, eq. (|25| ) is 
the generalization of the classical equation for the dispersion of surface waves in an elastic 
medium - the so called "Rayleigh waves" (see ref. ||12|| ) - to which it reduces when the 



dimensionless parameters X and Z, measuring the strength of the external electric and 
magnetic fields relative to the shear elastic forces, are equal to zero. 

The solutions of eq. ( P5| ) will be discussed in the next Section for various physical regimes. 
Before doing that, however, it is necessary to clarify a delicate point which arises when one 
attempts to take the limit C t — > of the above theory. Physically this corresponds to the 
very relevant case of a genuine liquid system, which is expected to have a vanishing shear 
modulus at low frequency. It is evident from eq. ( p9|) that this limit is singular: the variable 
£ tends to infinity, implying that only the u; component of the solution (that is, the part 
proportional to £ 2 in eq. (|29"D) survives. Nevertheless, the presence of a rapidly oscillating 
component of the solution with wavevector A ~ iq\£\, which does not vanish for x — > — oo, 
is disturbing. This difficulty becomes more evident when one tries to solve the equation for 
the dispersion: in the limit C t — > 0, £ — > oo, the argument of the square root lies on the 
negative real axis, where the presence of the branch cut prevents us from finding a solution. 

The resolution of these difficulties lies in the following physical considerations. Every 
system, liquid or solid, has a non vanishing shear viscosity r] at finite frequency. The shear 
viscosity contributes an additional term to the stress tensor (see [13||), which can be accom- 



modated within our formalism simply through the replacement /i — > fi — itur/ or, equivalently, 

C 2 t -> C 2 = C 2 - iuv, (30) 
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where v = rj/po is the so-called kinematic viscosity |T3[ . In a solid Ct remains finite for 



uj — > and therefore 1/ can be safely neglected. In a liquid, however, C t tends to faster 
than u, while v remains finite : the low frequency limit is therefore dominated by the 
viscosity. It is easy to see that the inclusion of viscosity eliminates the singularity of the 
solution in the limit Ct — > 0, v —>■ 0. This is because, after the replacement indicated in 
eq. (|30|), the wavevector qy/l — £ 2 in the limit C t — > reduces to ^J—iu/u = ^Ju/2u(l — i) 
which, for small u, has a very large positive real part, and therefore vanishes very rapidly 
away from the edge. Correspondingly, a solution of eq. ( |25| ) can always be found (at small 
q) if the viscosity is added according to the substitution indicated in eq. (|30D . In practice, 
the viscosity is expected to be small. In the following section we shall present our results 
treating u, C 2 , g, and uj c as formal parameters. The actual values of these parameters will 
be discussed in the concluding section. 



IV. THE DISPERSION RELATION IN VARIOUS PHYSICAL REGIMES 

The behavior of the solutions of eq. (|25| ) as functions of the dimensionless parameters X 
and Z is rather complicated. It is convenient to distinguish three different regimes according 
to whether the magnetic field, the edge electric field (g), or the shear modulus dominates. 

• I. Strong Magnetic Field 

This regime is characterized by \q\ <C \u c /C t \ and \q\ -C \oJ 2 /g\ Therefore \Z\ 3> 1 and 
\X\ <C \Z\ 2 . It is easy to see that eq. ( |25l) has two solutions in this limiting case. The 
first solution is 



x i x 2 Ji-xyz 2 -i 
W ^-z + z(z^- 4)L -xW 4| (31) 

The explicit form of the dispersion depends on the value of X. The most interesting 
case is \X\ 3> \Z\ which corresponds to long wavelengths |g| <C g/C 2 and vanish- 
ing shear modulus uo c C t g and is appropriate for the liquid state. In this case, 
substituting the definitions ( p6|) in eq. (^l|) one obtains, after simple manipulations 

gq ACfq\q\ g 2 \q\q g\q\^ 
uj = 1 ^ 4:ii/ — — , (32) 

UJ C UJ C 0J* U) A C 



where we have used eq. (|30| ) to obtain the imaginary part of the frequency. 

Thus, at sufficiently long wavelength, the dispersion is linear and independent of the 
visco-elastic constants. 

Recalling the definition of the "gravity acceleration" in eq. (fl|) we see that the phase 
velocity for long wavelength u/q coincides with the classical drift velocity v = cE/B, 
where E is the magnitude of the electric field at the edge. 

Other cases, compatible with our definition of the strong-field regime, can also be 
calculated from eq. (j3ip . For example, in the limit of vanishing electric field g — > or 
large shear modulus, such that |X| 1, \g/C 2 \ \q\ \uj c /C t \ we obtain 



S 



LU C UJ C U% 

Setting g = in this expression, we obtain the dispersion of surface elastic waves 
("Rayleigh waves") |L2j in the presence of a strong magnetic field. Notice that, due 
to the presence of the strong magnetic field, the dispersion is quadratic rather than 
linear in q. 

The second solution is 

(ii) i~Z + (X + A)/Z, (34) 

which implies 

gq 4:C^ 'q 2 

to = u c sgn(q) H 1 — sgn(q) — Aiuq 2 . (35) 

Notice that the two solutions obtained in this section describe waves propagating in 
opposite directions (i.e., the sign of the ratio uj/q is opposite in the two cases) and 
that the second solution, for q — > 0, is simply the manifestation of the uniform Kohn 



mode |TJ| at the edge of the system. 
II. Strong electric field (Gravity Waves) 

Let us now assume that the edge electric field is large, in the sense that io\j g <C |g| <C 
g/Cf. This implies \X\ ^> 1 and |X| 3> \Z\ 2 . Notice that q cannot tend to zero in this 
regime, unless u c = 0. The solutions of eq. are then found to be 



£~ ±\JX ±ZVX + A, (36) 



uj = ±^Jg\q\ ^1 + ^-^J + ^ - 2ivq 2 (37) 

It is comforting to notice that, for zero magnetic field and shear modulus, we recover 
the classical results for the dispersion and damping of "gravity waves" on the surface 
of a liquid [13| . 



III. Large Shear Modulus (Rayleigh Waves) 

Next, consider the case that the shear modulus dominates, in the sense that |g| ^> 
uJcjCt and |g| ^> g/C 2 . This implies |X| <C 1 and \Z\ <^ 1. Again, this definition does 
not allow one to take the limit q — > unless g and uj c are both zero. 

The doubly degenerate solution of eq. ( [25] ) in this case is £ = ±0.955, which implies 

(forg«yq>), 

u = ±0.955C t |g| - z(0.456)z/g 2 . (38) 



This is nothing but the familiar Rayleigh wave on the surface of an elastic solid [12 
with a small damping due to viscosity, and no correction from external electric and 
magnetic fields. 
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• IV. Zero Shear Modulus Limit (liquid state) 

In this case (and for zero viscosity) the equations of motion can be solved exactly. The 
expression for the two branches of the dispersion relation is: 




For q <u c 2 / g and g > w^/^ we recover respectively the expressions ( |3"2] ) and fl37|). 
The corresponding eigenfunction is given by u\ in eq. (|19|). 

We have thus exhausted all the physically different regimes. In case I to III, the form of 
the eigenfunction is obtained by substituting the appropriate value of £ in eq. (|29|) . Results 
for the lower branch of the complete numerical solution of eq. fl2"5|) exhibiting crossovers 
between the various regimes discussed in this section are shown in Fig. 1 and 2. In Fig. 1, 
the real part of the phase velocity \uuj c / ' gq\ is plotted against qC t /u) c for decreasing values 
of the ratio R = g/CtUJ c (R = 1,0.7,0.5) that is for increasing magnetic field or shear 
modulus or decreasing electric field. The Figure shows clearly the universal long wavelength 
behavior \ujuj c / gq\ « 1 for very small q. For R = 0.5 (strong magnetic field and finite shear 
modulus)and small q (q <C u) c /C t ) the exact curve agrees with the approximate expression 
eq. ( |3"3"D (labelled in the Fig. as "lw"). For R = 1 (strong magnetic field and small shear 
modulus) the third term in eq. ( |32"D becomes predominant and the exact curve lies below 1. 
For large q the exact curves approach in all cases the "Rayleigh Waves" regime, that is the 
lines \u>u c / gq\ = 0.955 /R. 

In Fig. 2 we plot the real part of \uiuj c /2q 2 C 2 \ versus qCf/g for \q\ <C \u c /C t \ for R = 0.01 
(solid line). This figure clearly shows a crossover from the linear dispersion controlled by 
the electric field to the quadratic dispersion controlled by the shear modulus in the case of 
strong magnetic field (see eq.(|33|)). The dashed line is the approximate curve for the "linear 
regime" (|g| <C \g/C 2 \, \uiui c /2q 2 C 2 \ = \g/2qC 2 \), while the dotted line corresponds to the 
"quadratic regime" (\g/C 2 \ |g| <C \u c /C t \, \uuj c /2q 2 C 2 \ = 1). 

In Fig. 3 we plot \u/u c \ against qg/ui 2 for the case of zero shear modulus (eq. fl39|)). 
For the low frequency branch (labelled with "— "), we plot also the leading term for long 
wavelength (\u/u c \ = qg/u 2 , dashed line) and the short wavelength approximate behavior 
(eq. (|3), dotted line). 

It is interesting to notice that no acceptable solution is found in the transition region 
between regimes II and III, that is, for uj c /Ct^q ~ g/Cf. 



V. DISCUSSION 



Up to this point our classification of different physical regimes has been purely formal: 
we have not yet specified the values of the parameters. We now wish to state the concrete 
predictions of our theory for typical systems at high magnetic field. 

The value of the electric acceleration at the edge is easily estimated as 

e 2 n d 

g = , 40 

m a 
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where no = po/m is the equilibrium density, d is the range of the interaction, and a is the 



width of the edge, which is typically of the order of the magnetic length I = ^hc/eB. In a 
magnetic field, it is convenient to introduce the filling factor u = 27rno/ 2 . Then our estimate 
for g takes the form 

eV d 
ZTch a 

The calculation of the shear viscosity is considerably more difficult. A mode-coupling 
calculation for the two-dimensional electron gas at zero magnetic field [|15j yields numerical 
results that can be accurately described by the formula [0 

v ~ (59r7 3/2 + cirj 1 + c 2 r7 2/3 + c 3 r" 1/3 ) _1 — , (42) 
v J m 

-,2 



where r s = Jl/itriQajQ is the usual electron gas parameter, a® = % /me 2 is the effective Bohr 
radius of the host semiconductor, and c = 0.25, C\ = 20.6, c 2 = 22.7, c 3 = 12.8. We are not 
aware of any calculation of the viscosity in the presence of a magnetic field, but we expect 
eq. fl42|) to give at least the right order of magnitude at a given density. 

As for the shear modulus, we expect it to vanish, if the system is liquid, at low frequency, 
leaving us with C 2 = —iujv. If, instead, the system is a solid, then C 2 has a finite value 
which can be estimated from dimensional considerations: 

2 V 2 

C? ^ (43) 
m 

It is the viscosity that can be neglected in this case. For example, in the classical two- 
dimensional Wigner crystal with Coulomb interactions (d = oo), the exact result is Cf 
0.2Ae 2 n /2 /m TL 



Let us first consider the dispersion of edge waves in a liquid. Because C t — 0, and v 

is very small we immediately see that \X\ ^> 1 and \Z\ ^> 1 at all realistic wavevectors 
q <C 1/a ~ 1/d. The inverse length uj 2 /g is given by 

^-«o'-?m 2 . (44) 



g v d\ I J 

This is of the order of \ ja or larger for typical quantum Hall systems in which the density 
is no ~ 10 10 — 10 11 cm -2 and the filling factor is less than 1. Therefore, these systems fall 
within the "strong magnetic field regime" of Section IV. Because the frequency vanishes 
linearly with q the ratio X/Z = gjuj c \J — iujv tends to infinity, and the dispersion is therefore 
given by eq. ( p2|) 

9Q , 9 2 q\q\ g\q\ 3 fA ^ 

to = 1 ^ Atu — — . (45) 

The eigenfunction has the simple form 

u oc [x + tsgn{q)y}e lqlx+iqy , (46) 
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which describes a circular motion of each volume element. Then, making use of our estimate 
(41) for the electric acceleration the phase velocity of the wave v = g/u c = (e 2 vo/h)(d/a) is 
obtained at once. 

Let us now consider the case that the low-frequency shear modulus does not vanish: 
this would happen, for example, if the electrons solidified in a Wigner crystal structure. 
According to our general discussion it might be possible, with increasing wavevector, to 
cross over to an "intermediate wavelength regime " , in which the dispersion is controlled by 
the shear modulus, and varies as q 2 (see eq. (p3|)). 

Unfortunately, this crossover is not likely to occur within the region of wavevectors in 
which our theory applies, namely q « 1/a ~ 1/d. Indeed, from eqs. (f4"0P , (f41~|), and ( f4~3|) 
we see that, for d comparable to, but somewhat larger than the average distance between 
electrons, both the ratios g/C 2 , and uj c /Ct are of the order of the inverse of the inter-electron 
distance. This implies that for q <C 1/a, 1/d both X and Z are 3> 1, and we clearly fall into 
the "long wavelength regime" of the previous section. The only possibility to observe the 
crossover to shear-modulus- sustained waves in the present model, would arise if the range 
of the interaction d were much less than the typical inter-electron distance - not an easily 
realizable situation. 

Thus, in summary, we have shown that the two-dimensional electron gas on a neutraliz- 
ing background of charge sustains only one macroscopically charged collective mode which 
decays exponentially as one moves away from the edge. The scale of this exponential decay 
is the same as the scale of variation of the density along the edge. We have derived analyti- 
cal and numerical expressions for the dispersion relations and the eigenf unctions in various 
physically distinct regimes. An important result of our study is that the character of the 
edge excitations in this model is controlled almost exclusively by the strength of the electric 
field at the edge, and does not depend significantly, on the shear modulus: therefore, the 
dispersion is the same, to leading order in q, for the liquid and the solid state (see eq. (|3TD ). 

These results suggest that the behavior of the I — V tunneling characteristics at low bias 
voltage would be the same at the edge of a liquid and of a Wigner crystal: in particular, the 
power law I ~ V 1 ^ , where uq is the bulk filling factor, is expected in both cases ||. 
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APPENDIX A: EXTENSION TO LONG-RANGE INTERACTION 

In this paper we have taken advantage of some mathematical simplifications arising from 
the assumption that the electron-electron interaction is of finite range in space. Extending 
the theory to properly include Coulomb interactions is nontrivial. In this Appendix we want 
to sketch an approximate method to do this extension, which entails minimal changes in the 



12 



structure of the equations. The method is non rigorous, yet it yields qualitatively correct 
results for the long wavelength dispersion. 

To begin, we observe that, in the equation of motion (|5|) the term C 2 V(V ■ u) should 
be replaced (ignoring the small contribution from K) by 



2 

V / p^— V' • u(r')dr'. (Al) 
J |r — r'| 

In an infinite system this term would pose no problem: by Fourier transformation it 
could be recast in the form of eq. ([5]) with a g-dependent longitudinal "sound velocity" 

= (A2) 
mq 

which tends to infinity as q tends to zero. In a semi-infinite system, however, things are not 
so simple, and the Coulomb interaction cannot be simply absorbed in a g-dependent sound 
velocity. 

The idea of our approximation scheme is to neglect the effect of the edge on the bulk 
equation of motion, which therefore retains the form of eq. (|5|), with C 2 tending to infinity 
in the long wavelength limit. The effect of the edge will be taken into account only via the 
boundary conditions ( [TT]) and (|12l), which force the solution to be a certain superposition 
of bulk waves. 

So far the theory is formally identical to the short-range case. However, observe that the 
key quantity g, which enters the boundary conditions, is ill-defined in the case of the Coulomb 
interaction. The external potential in this case is simply the electrostatic potential created 
by a uniform distribution of positive charge precisely compensating the electronic charge in 
the half-plane. The electric field produced by this charge tends to infinity (logarithmically) 
at the edge of the half-plane, and therefore g is infinite (see eq. ([[])). 

In order to obtain the correct form of the boundary conditions in the Coulomb case we 
return to eq. @ and note that from the equilibrium condition a x °J (r) is given by 

= -"3/ /, n uT^r', (A3) 

J Jx'<0 I — I | 

where the integral is restricted to the half-plane x' < 0. 

Consider a point in the vicinity of the geometric edge with r = (e + u x (0, y),y), where e 
is a length of the order of the physical width of the edge a. Similarly to what we did in the 
short range case in sec. II we now expand a x °] to first order in u x and discard the constant 
zero order contribution. After integrating with respect to x', we obtain 



oo ^2 



a x ° x \u x ,y) ~r%r & u x (0,y')dy'. (A4) 

J -°° ^e 2 + (y-y') 2 

Specializing to solutions with definite wavevector along the edge (i.e., u x (0, y) cx e iqy ) we 
see that eq. (|A4j) takes the simpler form 

*£>(u x ,y) ~ n 2 2e 2 K (\qe\)u x (0,y), (A5) 
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where Kq is the modified Bessel function. Unfortunately, the e — > limit of this expression 
does not exist, due to the logarithmic divergence of the Bessel function K (x) ~ —ln(x) for 
x — > 0. However, the divergence is very weak, and, in view of the fact that the position of 
the edge is defined within an uncertainty of the order of a (~ in our theory) it is legitimate 
(with logarithmic accuracy) to replace K (\qe\) by K (\qa\). 

Substitution of this expansion in eq. @ leads to the new boundary conditions: these 
can still be written in the form of eqs. (P%D, but now both the longitudinal sound velocity, 
and the electrical acceleration g are functions of q: the former is given by eq. (|^2|), while 
for the latter 

g^g(q) = 2e 2 ^K (\qa\). (A6) 
m 

From this point on, all the calculations proceed as in the short-range case. In particular, 
the logarithmic divergence of the effective g causes the dispersion of the edge waves to vary 
as — (2z/ e 2 //i)gln(ga) in a magnetic field, and as y , 2e 2 n ^ 'qln(qa) without a magnetic field. 

It is amusing to observe that in the three-dimensional case the effective g at the surface 
would be g(q) = 2iTpoe 2 /m 2 q , leading, in the absence of a magnetic field, to the well known 
result for the frequency of surface plasmons, namely u = ^fgq = J2iie 2 n /m. 
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FIGURES 



FIG. 1. Real part of the dimensionless phase velocity \uu c /gq\ versus qCt/co c for three values 
of R = g/CtuJc- For each curve the short wavelength limit \uu c /gq\ = 0.955/i? - coinciding with 
the dispersion of classical Rayleigh waves - is plotted for comparison. For the curve with R = 0.5 
we also plot the long wavelength limit \uuj c /gq\ = 1 + 2C 2 \q\/g ("lw" curve). 

FIG. 2. Real part of \uJUJ c /2q 2 C 2 \ versus qC 2 /g for R = 0.01 (solid line). The dashed line is 
the approximate curve for the "linear regime" (\q\ <C \g/Cf\, \iouj c /2q 2 C 2 \ = \g/2qCf\), while the 
dotted line corresponds to the "quadratic regime" {\g/C 2 \ <C \q\ <C \u c /Ct\, \couj c /2q 2 C 2 \ = 1). 

FIG. 3. Real part of the upper ("+") and lower ("— ") branch of the dimensionless dis- 
persion curve \u/u c \ versus qg/io 2 for vanishing shear modulus. The dashed line represents 
the long wavelength limit (|u;/a; c | = qg/m 2 ) and the dotted line the short wavelength limit 
(\u}/u c \ = v 'qg/^'c ~ 1/2) of the lower branch. 
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